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We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling 
limits from the recent results of the LUX direct search dark matter experiment. After a brief review 
of WIMP-nucleon scattering, we derive LUX limits on the spin-dependent WIMP-nucleon couplings 
over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical 
parameters. We find that, under these and other common assumptions, LUX excludes the entire 
spin-dependent parameter space consistent with a dark matter interpretation of DAMA’s anomalous 
signal, the first time a single experiment has been able to do so. We also revisit the case of spin- 
independent couplings, and demonstrate good agreement between our results and the published 
LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory 
in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via 
axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to 
place constraints on generic dark matter theories. 


I. INTRODUCTION 

Evidence for a large amount of non-baryonic “dark 
matter” (DM) in the universe has been accumulating for 
decades mm- Recent observations of the cosmic mi¬ 
crowave background have provided a precise measure¬ 
ment of the DM relic density, and also strongly sup¬ 
port the idea of “cold” DM in the form of weakly in¬ 
teracting massive particles (WIMPs) [3H5]- The failure 
of the Standard Model (SM) of particle physics to ad¬ 
equately explain a variety of astrophysical observations 
has prompted the development of a large number of par¬ 
ticle theories beyond the SM [5]. 

Concurrently with these theoretical developments, a 
large number of experiments have been conducted to 
search for DM annihilation in distant astrophysical ob¬ 
jects, produce and observe DM particles in high energy 
particle collisions, or observe the direct interaction of par¬ 
ticles of DM with Earth bound detectors. Although there 
are tantalising hints of DM signatures in one or more of 
these experiments, there is as yet no uncontroversial de¬ 
tection of (non-gravitational) WIMP interactions with 
ordinary matter [7]. Given a particular new theory of 
particle physics, it is sometimes challenging to assess the 
likelihood of the model (as a function of the model pa¬ 
rameters) given the null results of these experiments. 

In this paper, we present LUXCalc, a new utility for 
assessing the likelihood of new physics models given the 
recent results of the LUX experiment [8], the most con¬ 
straining direct search experiment to date for a wide 


range of WIMP models. In addition to describing the 
use of the tool for the general user, we apply it to derive 
LUX limits on spin-independent (SI) and spin-dependent 
(SD) WIMP-nucleon couplings. The former show good 
agreement with the official LUX results, thus validating 
our approach. Whilst SD limits have not been provided 
by the LUX collaboration, our results in the neutron-only 
coupling case are in close agreement with those calcu¬ 
lated by Ref. [5] We present for the first time the LUX 
SD proton-only limits and discuss the general SD mixed 
coupling case, finding that LUX alone fully excludes the 
SD-interacting DM interpretation of the anomalous sig¬ 
nal seen in DAMA m- 

The paper is structured as follows. In Section [IT] we 
provide a brief and self-contained review on the physics of 
WIMP-nuclear scattering, in both the SI and SD case. In 
Section m we review the LUX experiment, and present 
the methodology for our determination of the LUX likeli¬ 
hood and constraints. Section [rv| presents our SI and SD 
coupling limits as a function of the WIMP mass, includ¬ 
ing comparisons with the published LUX limits (where 
relevant), and those of other DM experiments. We also 
place constraints on an effective DM theory for which 
the SD constraints are particularly important: the case of 
fermionic DM interacting with a spin-1 mediator that has 
an axial-vector coupling to both the WIMP and to SM 
fermions. Pinally, we present conclusions in Section [V] 
The Appendix describes the use of the LUXCalc software. 
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II. THEORY REVIEW 


Direct detection experiments aim to observe the recoil 
of a nucleus in a collision with a DM particle [H]. After 
an elastic collision with a WIMP x of mass , a nucleus 
of mass M recoils with energy E = (/r^u^/M)(l — cosd), 
where /i = m^M/{rriy^ + M) is the reduced mass of the 
WIMP-nucleus system, v is the speed of the WIMP rel¬ 
ative to the nucleus, and 9 is the scattering angle in the 
center of mass frame. The differential recoil rate per unit 
detector mass is 


dR uq / da \ 
dE ~ M VdE/ 



da 

dq^ 


(I) 


where Uq = pojm^ is the number density of WIMPs, 
with po the local DM mass density; /(v,t) is the time- 
dependent WIMP velocity distribution; and 
is the velocity-dependent differential cross-section, with 
= 2ME the momentum exchange in the scatter. In 
the typical case that the target material contains more 
than one isotope, the differential rate is given by a mass- 
fraction weighted sum over contributions from the iso¬ 
topes, each of the form given by Eqn. Q. Using the 
form of the differential cross-section for the most com¬ 
monly assumed couplings. 


dR 


I 


o-(9) min 


( 2 ) 


where a{q) is an effective scattering cross-section and 


, j3 /(vU) 

r][Vmin,t)= / dv - (3) 

is the mean inverse speed, with 


Vmin = 



(4) 


the minimum WIMP velocity that can result in a re¬ 
coil energy E. Equation Q conveniently factorizes the 
differential rate into particle physics terms {a) and astro¬ 
physics terms {pop), which we describe separately in the 
following sections. More comprehensive reviews of direct 
detection can be found in Refs. [aiiMi]- 


A. Particle physics: cross-section 

For a SUSY neutralino and many other WIMP can¬ 
didates, the dominant WIMP-quark couplings in direct 
detection experiments are the scalar and axial-vector cou¬ 
plings, which respectively give rise to SI and SD cross- 
sections [11]. In both cases, 

( 5 ) 


to leading order. Here, 0 is the Heaviside step func¬ 
tion, qmax = 2/ru is the maximum momentum transfer in 
a collision at a relative velocity v, and the requirement 
q < 9max gives rise to the lower limit v > Vmin in the in¬ 
tegral for T/ in Eqn. (|^. In the above equation, ag is the 
scattering cross-section in the zero-momentum-transfer 
limit (we shall use crsi nnd crgD to represent this term in 
the SI and SD cases, respectively) and E^(q) is a form 
factor to account for the finite size of the nucleus. The 
WIMP coherently scatters off the entire nucleus when the 
momentum transfer is small, giving F^{q) —>■ 1. However, 
as the de Broglie wavelength of the momentum transfer 
becomes comparable to the size of the nucleus, the WIMP 
becomes sensitive to the spatial structure of the nucleus 
and B^(g) < I, with B^(g) <C I at higher momentum 
transfers. It is traditional to define a form-factor cor¬ 
rected cross-section 


a{q) = agE'^iq) , (6) 

as was used in Eqn. (§ above. We note that this is 
an effective cross-section, whereas the actual scattering 
cross-section is given by / dq'^-^{q^,v) for a given rela¬ 
tive velocity v. The total WIMP-nucleus scattering rate 
is then the sum over both the SI and SD contributions, 
each with its own form factor. 


1. Spin-independent cross-section (SI) 


The SI WIMP-nucleus interaction, which occurs 
through operators such as {xx){dd)^ lias the cross-section 


CTSI = 


ZGIMA-Z)GI, 


4^2 r i2 

= Zff+{A-Z)U 

TT L J 

(7) 

where Z and A — Z are the number of protons and neu¬ 
trons in the nucleus, respectively, and /p (/„) is the effec¬ 
tive coupling to the proton (neutron), with the alternate 
normalization = 2ff^ also found in the literaturej^ 
For neutralinos and most other WIMP candidates with 
a SI interaction arising through scalar couplings, /„ ~ /p 
and the SI scattering cross-section of WIMPs with pro¬ 
tons and neutrons are roughly comparable, (Tn,si « o'p,Si- 
For identical couplings (/n = /p), the SI cross-section can 
be written as 


csi = cTp SI, 


( 8 ) 


where /ip is the WIMP-proton reduced mass. As neu¬ 
tralinos are the currently favored WIMP candidate, this 
assumption is widely made throughout the direct detec¬ 
tion literature, though models can be constructed that 


2 Notably, are the Gir-like effective four-fermion coupling con- 
stants in the case of scalar interactions m and are the normal¬ 
ization used in DarkSUSY |17| . MicrOMEGAs uses Ajy = § | 18| . 












3 


violate this /„ — /p condition (e.g. isospin-violating DM 
m)- We will consider only identical SI couplings case 
when later examining the LUX results. 

The SI cross-section grows r^idly with nuclear mass 
due to the A'^ factor in Eqn. which arises from the 
fact that the total SI coupling of the WIMP to a nucleus 
is a coherent sum over the contributions from individual 
protons and neutrons within. Direct detection experi¬ 
ments therefore often use heavy nuclei to increase their 
sensitivity to WIMP scattering. 

The SI form factor is essentially a Fourier transform 
of the mass distribution of the nucleus. A reasonably 
accurate approximation is the Helm form factor naEDi; 

F(,) = cos(,r„) _ 

\Vn) 

where s ~ 0.9 fm and -I- — 5s^ is an effective 

nuclear radius with a ~ 0.52 fm and c ~ 1.23A^/^ — 
0.60 fm. Further details on SI form factors can be found 
in Refs. [nm]. 


S. Spin-dependent cross-section (SD) 


The SD form factor is given in terms of the structure 
function S'(g) as 


F^{q) 


S{q) 

S{Q) 


( 11 ) 


such that F^(0) = I as previously prescribed. The S{q) 
have the functional form 


S{q) = Op S'pp (g) + a2 5n„(g) + apa„^p„(g) . (12) 

The differential SD scattering cross-section can alterna¬ 
tively be written in terms of this structure function a^ 


In the limit of zero momentum transfer {q —>■ 0) the func¬ 
tions SxyiO) {x,y = p,n) can be related to expectation 
values of the proton/neutron spin [2^1^ : 


(^p)^ 


(J+1) (2J+1)^PP^°^ 
(J+1) (2J+1)^“^°^ 


and 


(14) 


SD scattering is due to the interaction of a WIMP 
with the spin of the nucleus through operators such as 
(X7 ai75X)( 97 ^ 759 )I takes place only in those detec¬ 

tor isotopes with an unpaired proton and/or unpaired 
neutron. The SD WIMP-nucleus cross-section is 


ctsd = 


4/r2 (j+l) 


TT J 
32^^ (J-bl) 

TT J 


GL(*5p) + G",(5„) 


Up(aSp) -f Un(‘Sn) 


( 10 ) 


where Gp is the Fermi constant, J is the spin of the nu¬ 
cleus, (Sp) and (Sn) are the average spin contributions 
from the proton and neutron groups, respectively, and 
Op (ttn) are the effective couplings to the proton (neu¬ 
tron) in units of 2y/2GF; the alternative normalization 
Gf, = 2V2Gp apf is also found in the literature]^ Un¬ 
like the SI case, the two SD couplings Up and Un may 
differ substantially (though they are often of similar or¬ 
der of magnitude), so that a simplification comparable 
to Eqn. (j^ is not made in the SD case. Because of 
the uncertain theoretical relation between the two cou¬ 
plings and following from the fact that one of {Sp) or 
(S'n) is often much smaller than the other, experiments 
typically only significantly constrain one of the two SD 
cross-sections, Up^D or Cn.SDj but not both. 


Under an alternative basis 


ao = Up -I- On Oi = Op - On , (15) 

which is a more convenient basis for nuclear physics work, 

S{q) = alSoo{q) -{-aj Sii{q) -f oooi S'oi(g) , (16) 

where the structure functions in the two bases are related 
by 


Spp = Sqq + Sii -\- S'oi 

I^nn = Soo + S'n — Soi (17) 

Spn = 2 (Soo ~ Sii) . 

The Sij {i,j = 0,1) are calculated in the literature. For 
the two non-zero-spin xenon isotopes (^^®Xe and ^^^Xe), 
we take these structure functions from Ref. [21] , utilizing 
one plus two body (lb-f2b) axial-vector currents where 
possible. 


B. Astrophysics: dark matter distribution 

The DM halo in the local neighborhood is most likely 
dominated by a smooth and well-mixed (virialized) com¬ 
ponent with an average density po. The simplest model 


® Though we use ajv and to distinguish between the two nor¬ 
malizations here, ajv is frequently used within the literature for 
both cases. The G^ normalization is used by DarkSUSY dsiin], 
while micrQMEGAs uses = US]- 


The factor of is omitted if Eqn. is written in terms of 
rather than a^v- Again, sometimes the latter notation is 
used in the literature to represent the former normalization as 
defined here. 
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for this smooth component is often taken to be the Stan¬ 
dard Halo Model (SHM) [25l[26], a non-rotating isother¬ 
mal sphere with an isotropic, Maxwellian velocity distri¬ 
bution and most probable speed vq, where for the SHM, 
vq is equal to the disk rotation speed Vrot- The SHM ve¬ 
locity distribution in the Galactic (i.e. non-rotating) rest 


frame is 




/(v) = 1 

\ 

[0. 

\3/2 2 2 

j for |v| < iiesc 

otherwise. 

(18) 

Here, 

N — 

■^’esc — 

2 2 

erf( 2 :) - ^ze ^ , 

y/TT 

(19) 

with 0 = 

^esc/ 

is a normalization factor. 

The 


Maxwellian distribution is truncated at the escape ve¬ 
locity Vesc to account for the fact that WIMPs with suffi¬ 
ciently high velocities escape the Galaxy’s potential well. 

Ultimately, the velocity distribution in the Earth’s rest 
frame is the one relevant for direct detection, which can 
be obtained after a Galilean boost: 

/(v) = /(v© -b v), (20) 

where v© = vlsr + v©_pec is the motion of the Sun rela¬ 
tive to the Galactic rest frame, vlsr = (0,Vrot,0) is the 
motion of the Local Standard of Rest in Galactic coor¬ 
dinates!^ and v©_pec is the Sun’s peculiar velocity. The 
additional time-dependent orbital motion of the Earth 
about the Sun is postulated to give rise to a measurable 
modulation in the signal m [26] : indeed, some experi¬ 
ments have claimed positive results for such signatures 
m- This small Earth orbital motion and resulting time 
dependence has been neglected here as it is not relevant 
to the LUX calculations. 

The choice of values for the various DM distribution 
parameters is important for interpreting direct detection 
results. Due to the inability to directly observe the DM 
and various systematic issues in interpreting observations 
of the Galaxy, some of those parameters are known to 
limited precision. The range of estimates for these pa¬ 
rameters are shown in Table |l] as well as values commonly 
used for comparing direct detection results (“canonical”) 
and the default values used by LUXCalc. The LUXCalc 
default values for po and Vrot are somewhat larger than 
the historical canonical values as more recent estimates 
tend to prefer the somewhat larger values. However, the 
canonical values are not inconsistent with the current 
observations and will thus continue to see wide usage. 
The parameter values can be easily changed in LUXCalc 
and we will use the canonical values ourselves when com¬ 
paring LUX results with other experiments in later sec¬ 
tions. In addition to these SHM parameters, we take 
v©,pec = (ll,12,7) km/s[39]. 


® Galactic coordinates are aligned such that x is the direction to 
the Galactic center, y is the direction of the local disk rotation, 
and z is orthogonal to the plane of the disk. 


III. THE LUX EXPERIMENT AND ANALYSIS 

In this section, we describe the basic operation of 
the LUX DM search detector and how their data is 
used to constrain WIMP parameters. The LUX ex¬ 
periment uses a liquid-xenon time projection chamber 
(TPC) to identify DM recoil events and distinguish them 
from other (background) events [1D|. The first LUX re¬ 
sults, released in 2013 [S], involved a fiducial exposure 
of 1.01 X 10^ kg-days, comparable to that of the then- 
leadin^XENONlOO experiment [41] (7.6 x 10^ kg-days), 
which likewise used a TPG detector with a xenon target. 
The slightly larger exposure and somewhat better detec¬ 
tor performance characteristics allowed LUX to overtake 
XENONIOO and LUX is now the leading experiment in 
terms of sensitivity to a large variety of WIMPs (by about 
a factor of two over XENONIOO). 

The principles of operation for a TPC are as follows. 
A recoiling xenon nucleus (or recoiling electron, in the 
case of some background processes) in the liquid xenon 
target will interact with other nearby atoms, inducing 
both excitations and ionization of those atoms. The rela¬ 
tively quick relaxation of the excitations releases photons 
that are collected by photomultiplier tubes (PMTs); this 
prompt scintillation is referred to as the SI signal. While 
some of the ionized electrons can recombine (possibly 
producing excitations and contributing to the SI prompt 
scintillation signal), many of the free electrons are drawn 
away from the interaction site to the surface of the liquid 
by the application of an electric field across the liquid. 
Above the liquid is a small region of gaseous xenon un¬ 
der a higher electric field. Electrons reaching this region 
rapidly accelerate and collide with xenon atoms in the 
gas, releasing photons and more electrons in a cascade 
process; this “proportional scintillation” is the S2 signal. 
The drift time of the ionized electrons across the liquid 
is substantially longer than the relaxation time of the 
xenon excitations, so the SI and S2 scintillation signals 
are easily distinguishable. The benefit of observing two 
signals is that electron recoils, produced by background 
radiation, tend to produce a relatively larger amount of 
S2 for a given SI, so S2/S1 is used as a background dis¬ 
crimination parameter. 

For a given WIMP spectrum, the average expected 
number of signal events in some analysis region is 

rOO 1JD 

= dE (t>{E)—{E), (21) 

where MT is the detector mass x time exposure and ^(E) 
is the fraction of recoil events of energy E that will be 
both observed and fall into the predefined analysis re¬ 
gion. This ((>(i?) detection efficiency accounts for various 
trigger efficiencies, intrinsic statistical fluctuations and 


In terms of sensitivity to WIMPs with SI interactions and masses 
above ~10 GeV. 
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estimate 

canonical 

LUXCalc default 

Refs. 

po [GeV/cm®] 

0.2 - 0.7 

0.3 

0.4 

[28H33] 

Wrot [km/s] 

200 - 250 

220 

235 

[SUIT] 

Veac [km/s] 

500 - 600 

544 

550 

m 


TABLE I. Values of parameters describing the Standard Halo Model (SHM). The columns show in order: ranges of values 
typically found in the literature, commonly used values for comparison of experimental results, and the default values used by 
LUXCalc. In the SHM, the final parameter vq is not independent, but Hxed by the relation vo = Urot- 


the PMT response (i.e. detector resolution), and analysis 
cuts. The benefit of the above form is that all the com¬ 
plicated detector physics and responses are rolled into 
4>{E), which is independent of the type of WIMP inter¬ 
action or spectrum. Thus, for a given experimental result 
and analysis region, 4>{E) can be tabulated once and then 
used for analyzing arbitrary WIMPs. We use the TPCMC 
monte carlo [42] to model the detector response and gen¬ 
erate the relevant 4>{E) efficiency curves. TPCMC relies on 
NEST for modeling the microphysics of a recoiling 

xenon atom. 

We take as our analysis region in the S2/S1-S1 plane 
the region above the S2 > 200 PE threshold, below the 
nuclear recoil calibration data mean S2/S1 curve, and 
2 PE < SI < 30 PE. This region matches that used by the 
LUX collaboration, except for the imposition of the hard 
S2/S1 cut that is not necessary for their profile likelihood 
analysis [15]. This region contains one observed event at 
an SI of 3.1 PE with a mean expected background of 5 = 
0.64 events. Ideally, a lower S2/S1 bound to the region 
should be imposed in the count-based analyses we will be 
using, but such a bound can often be set low enough that, 
in practice, the (j){E) are negligibly affected. The lower 
bound would serve more to exclude very low S2/S1 events 
that are almost certainly backgrounds; luckily, there are 
no such events in the LUX results and this issue is mootQ 

Given an observed number of events N and expected 
background b — being 1 and 0.64 for this LUX analysis, 
respectively — a likelihood can be constructed from the 
Poisson distribution, with 

Cim^,a\N)=P{N\m^,cr)=^ ^ ( 22 ) 

where /r = ^(m^,cr) is the expected number of signal 
events for a given WIMP mass and one or more 
scattering cross-sections cr (or, alternatively, couplings). 
This likelihood can be easily combined with those from 
a variety of other physics data, such as from colliders or 
indirect DM searches, in statistical scans of the Minimal 
Supersymmetric Standard Model (MSSM) or other new 
physics frameworks; see e.g. [471151] . 


^ A conservative imposition of a lower S2/S1 cut corresponding to 
the 10% lower tail of the nuclear recoil calibration data (dashed 
red curve in Figure 4 of Ref. 0 ) results in constraints that are 
weaker by ~ 20-30%. 


We derive constraints in a cr-m^ parameter space 
via two methods: the Feldman-Cousins (Poisson-based) 
method |52] and Yellin’s maximum gap method [53] . 
These two techniques correspond to analyses with and 
without background subtraction, respectively, with the 
latter case commonly found in the literature due to the 
difficulty in reliably identifying and characterizing back¬ 
ground contributions. 

The Feldman-Cousins method derives a confidence in¬ 
terval in a for each m-^ (resulting in a raster scan in the 
cr-m^ plane) that is consistent with the observed number 
of events given the expected background. This confi¬ 
dence interval may be either one or two sided and, thus, 
is capable of excluding the zero-signal case when excess 
events are found. This method is relatively straightfor¬ 
ward and easy to implement, but one of the drawbacks is 
the lack of spectral information in the analysis: only the 
total counts are used, not the distribution of SI values 
(a coarse proxy for recoil energy). This spectral informa¬ 
tion can be useful in distinguishing between signal and 
background and can aid in constraining the mass of the 
WIMP in the event of a positive result, as heavier WIMPs 
induce more energetic xenon recoils and larger SI values, 
on average. Analyses that make use of spectral informa¬ 
tion are substantially more complex to implement and 
are thus typically only performed by the experimental 
collaborations themselves; see e.g. Refs. [54l[55] . 

The maximum gap method makes no presumptions 
about the amount of background that might be con¬ 
tributing to the observed events, instead assuming any 
or all of the events might be signal to generate a conser¬ 
vative exclusion limit in a at each m^: any cross-sections 
above this limit would yield too many events even if back¬ 
ground contributions were ignored. This method does 
make use of the SI distribution of the observed events, 
focusing on an SI range where the expected number 
of events is largest relative to the number of observed 
events. Specifically, the maximum gap method breaks 
the full observable (SI) range into intervals separated by 
the observed events. If /is, are the predicted number of 
events in each of these intervals, the “maximum gap” is 
the one where /ifc is largest. If x = is the fraction 

of signal events expected in this interval, then 


Co{x,fi) 


E 


k^O 


{kx- 



— kx J 


(23) 


kl 
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is the probability of the maximum gap being smaller 
than observed. In other words, a WIMP mass and cross- 
section(s) is excluded at greater than a 90% confidence 
level (CL) if Co > 0.9. To perform this calculation with 
the LUX result (with one event at an SI of 3.1 PE), we di¬ 
vide the previously discussed LUX S2/S1-S1 analysis re¬ 
gion into two parts, SI G [2,3.1] PE and SI G [3.1,30] PE, 
and use TPCMC to generate efficiency factors (j)i{E) and 
02 (E), respectively, for the two intervals]^ The expected 
number of events in an interval Hk can then be calculated 
via Eqn. (211 under the replacement 0 —>■ 0^. 

In tandem with this paper, we provide LUXCalc, a soft¬ 
ware package for performing the above-described LUX 
analyses. LUXCalc can be used as a standalone program 
or as a library to be called from other software packages 
such as DarkSUSY [T3 [SB] and micrOMEGAs A 

description of LUXCalc and its usage can be found in the 
Appendix. The software can be found at Ref. [SD] or as 
ancillary files to the arXiv version of this paper. 

Before turning to our results, we take a moment to 
stress that the 4>{E) and (j)kiE) curves used here require 
a full statistical modeling of LUX’s TPC detector to gen¬ 
erate and cannot be trivially generated from any efficien¬ 
cies provided by LUX in Ref. [8] . One might be tempted 
to take the no-S2/Sl cut efficiency from Figure 1 of that 
reference and apply an additional factor of 0.5 to account 
for the fraction of nuclear recoils falling below the S2/S1 
mean in the calibration data. There are two reasons why 
this is incorrect. First, the S2/S1 cut is not independent 
of the other cuts and, in fact, is very highly correlated 
with the SI G [2,30] PE cut near the boundaries. Second, 
the mean of the calibration recoil band in S2/S1-S1 space 
represents the distribution convolved over energy for that 
particular calibration spectrum; this does not imply that 
50% of the events at any specific energy will fall be¬ 
low that mean. Finally, the efficiency curve provided by 
LUX (without the S2/SI cut), applies only to the whole 
[2,30] PE SI interval and cannot be decomposed into the 
subintervals used for the maximum gap analysis. As we 
show in the next section, our efficiency curves nearly ex¬ 
actly reproduce LUX’s own low-mass constraints, an in¬ 
dication that those efficiency curves are correctly mod¬ 
eled (the naive approach of simply applying an additional 
50% nuclear-recoil-median-cut acceptance would lead to 
constraints too weak by a factor of two at low masses). 
A tabulation of these efficiencies is also included as an 
ancillary file on the arXiv version of this paper. 


IV. PHYSICS RESULTS 

In this section, we apply the previously described 
methods to analyze various physics models. We first ex¬ 
amine in Section |IV A| the LUX constraints on generic 



FIG. 1. Spin-independent cross-section constraints for the 
LUX experiment. Constraints are generated with and with¬ 
out background subtraction: a Poisson-based analysis with 
one observed event and 0.64 expected background events is 
used in the former case (dashed black), while the maximum 
gap method is used in the latter case (solid black); see the 
text for details. For comparison, the official LUX collabora¬ 
tion constraints are also shown (red), based upon an event- 
likelihood analysis. All constraints are at the 90% CL; cross- 
sections above the curves are excluded at greater than this 
level. 


WIMPs with SI and SD interactions, comparing our re¬ 
sults with those from other DM searches and from the 
LUX collaboration itself. In Section |IVB[ we then apply 
LUX constraints to an effective theory model where SD 
interactions are particularly relevant. 


A. Generic coupling limits 

Our LUX SI scattering constraints are shown in Fig¬ 
ure with the maximum gap limit shown in solid black 
and the Poisson-based constraint shown in dashed black; 
constraints are at the 90% CL. The official LUX limits 
are also shown (thin red). Though they are not our pre¬ 
ferred parameters, we use vg = Vrot = 220 km/s and 
Po — 0.3 GeV/cm^ to generate constraints in this and 
later figures to allow for a proper comparison with vari¬ 
ous other experimental constraints that use these values, 
such as the official LUX limits. The maximum gap limit 
is remarkably close to the LUX collaboration’s own limit, 
while the likelihood-based constraint is weaker by ~30% 
except at lower WIMP masses (which will be discussed 
momentarily)!^ There are two main potential reasons 


® While the likelihood analysis can produced two-sided limits, in 
this case the lower limit is simply zero cross-section, i.e. no 
WIMP signal: the LUX result is consistent with backgrounds 
alone at the 90% CL in this analysis. 


By definition, 0(£) = 0fe(E). 
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FIG. 2. Spin-independent cross-section constraints from 
various direct searches. Parameters that can reproduce the 
CoGeNT (90% CL), CRESST (2cr), and DAM A (3a) anoma¬ 
lous signals are shown In the filled regions. SuperCDMS and 
LUX exclusion constraints are shown at the 90% CL. The 
LUX maximum gap constraints (solid black) are shown with 
nuclear recoil spectra limited to E > 0,1,2,3 keV (strongest 
to weakest). The official LUX limit used a conservative 3 keV 
minimum. 


why the likelihood-based constraint is weaker: (1) the 
analysis region for our analysis is the lower-half of the 
nuclear recoil band and thus contains only about half of 
the potential signal (while the LUX collaboration analy¬ 
sis uses all of it) and (2) our Poisson likelihood makes no 
use of spectral information (i.e. event SI). The first rea¬ 
son is not, in fact, a major issue in practice as the upper 
half of the nuclear recoil band that is being ignored is 
contaminated by background events and does not signif¬ 
icantly improve the signal-to-noise in those analyses that 
use it. 

The LUX SI scattering constraints for low WIMP 
masses are shown in Figure For comparison, the 
parameters consistent with a DM interpretation of the 
anomalous signals seen in CoGeNT [61], CRESST [62] . 
and DAMA [TT] |B3] are shown, as well as the exclusion 
constraints from SuperCDMS [53] . As the WIMP mass is 
lowered, the Poisson-based likelihood constraint (dashed 
black curve) becomes comparable to and then slightly 
better than the maximum gap limit (farthest left solid 
black curve, other curves discussed below). The low-mass 
improvement of the likelihood case relative to the maxi¬ 
mum gap case is due to the fact that the single event in 
the analysis becomes consistent with the expected WIMP 
spectrum, leading to a slight weakening in the sensitivity 
of the maximum gap method. 

As measurements of the scintillation and ionization at 
very low nuclear recoil energies are limited and the theo¬ 
retical models will eventually break down at sufficiently 
low energies (see e.g. Ref. [5^), the LUX collaboration 


FIG. 3. Spin-dependent WIMP-neutron cross-section con¬ 
straints for the neutron-only coupling case. 


conservatively ignores contributions from WIMP scatter¬ 
ing events with E < 3 keV when analyzing their results. 
While this has little impact on constraints for WIMPs 
heavier than ~20 GeV, it becomes important for light 
WIMPs as light WIMPs can only induce low-energy re¬ 
coils. In the figure, we show our own maximum gap con¬ 
straints when considering only U > 0, 1, 2, and 3 keV 
(solid black curves from left to right, or most to least 
constraining). With the same E > 3 keV that LUX 
uses, the maximum gap constraint closely matches their 
constraint. To be clear, placing a minimum on the con¬ 
tributing E is not quite the same as defining the thresh¬ 
old in the detector. The actual trigger and SI analysis 
thresholds are already built into the (j>{E) efficiency term 
in Eqn. (211 regardless of the choice of lower bound in 
the integral over E. That efficiency falls rapidly below 
3 keV, from 23% at 3 keV, to 4% at 2 keV and 0.04% 
at 1 keV. Even if the lower bound of integration is set to 
E = 0 keV, no recoil events with energies below 0.5 keV 
will contribute to the signal as (j>{E) is zero at these en¬ 
ergies. Placing a minimum requirement on E serves to 
avoid the (already suppressed) contributions from events 
where the NEST model has little experimental data to 
ensure its accuracy. For our remaining analyses, we do 
not apply this artificial cut on the low-energy recoil spec¬ 
trum, though this choice has little impact on our conclu¬ 
sions and constraints with a cut applied can be easily 
generated with the LUXCalc code as described in the Ap¬ 
pendix. 

We now turn to the SD case, with WIMP-nucleon 
cross-section constraints shown in Figures [^ & [^ for 
neutron-only (op = 0) and proton-only (on = 0) cou¬ 
plings, respectively. As xenon has neutron-odd isotopes, 
LUX should be particularly sensitive to a WIMP with 
neutron-only SD couplings. In Figure]^ we show by the 
black curve the LUX constraints in this case, as deter- 










FIG. 4. Spin-dependent WIMP-proton cross-section con¬ 
straints for the proton-only coupling case, including indirect 
search limits from the IceCube/DeepCore experiment. Note 
the IceCube/DeepCore W^W~ constraint uses the tau chan¬ 
nel for masses smaller than that of the W boson. 


mined via the method described in Section Ell We also 
show constraints from other experiments with neutron- 
odd target materials: CDMS II [66] (at lower masses, 
also for the low-threshold analysis [57]), ZEPLIN-III [55] . 
and XENONIOO [55]. Here, the somewhat improved ex¬ 
posure and threshold of LUX over that of XENONIOO 
is evident by the ~ x2 stronger constraint at heavy 
WIMP masses and the extension of the constraints to 
lower WIMP masses before losing sensitivity. 

A SD proton-only interaction has historically been a 
means of producing DAMA’s anomalous signal while 
evading the null results of other searches [70] : the proton- 
even target isotopes used by many experiments have sup¬ 
pressed interactions (and thus little expected signal) in 
this case, while the proton-odd sodium-iodide target used 
by DAMA ensures they remain sensitive to the WIMP. 
Recent results from the SIMPLE [7T], PICASSO [75] . 
and COUPP [^ experiments, which also have proton- 
odd target isotopes, are now in conflict with a SD proton- 
only coupling explanation for the DAMA signal mm 
as shown in Pigure [^ Due to the ever-increasing de¬ 
tector sizes, even experiments with proton-even targets 
like XENONIOO [69] are starting to probe DAMA’s pre¬ 
ferred parameter region. Our determination of the LUX 
constraints are shown by the thick black curve. We see 
here that LUX, even though it uses a proton-even xenon 
target material, fully excludes the DAMA region. 

Indirect DM searches via neutrinos produced when 
WIMPs are caught and then annihilate in the Sun can 
place constraints on the SD WIMP-proton cross-section 
as collisions of WIMPs with hydrogen (protons) is part 
of the process for capturing WIMPs in the Sun m- Fig¬ 
ure]^ shows constraints placed by the IceCube/DeepCore 
m and Baksan m neutrino detector searches for such 


neutrinos. The constraints depend on the annihilation 
channel and are shown here for the representative 6-quark 
and W-boson channels. While neutrino searches can be 
very sensitive to WIMPs with SD proton couplings, the 
high thresholds in IceCube/DeepCore and some other 
neutrino experiments means they are often unable to 
probe for light WIMPs as LUX and other direct searches 
are capable of doing. Furthermore, the limits shown here 
assume the DM capture and annihilation processes in the 
Sun are in equilibrium, an assumption that may not be 
true for many DM candidates [77]. 

The exclusion of the DAMA region by LUX in the SD 
proton-only coupling case, the case where LUX limits are 
approximately at their weakest, suggests that the LUX 
result may exclude any SD explanation for the DAMA 
signal, the first time a single experiment would be able to 
do so. After a more careful examination over the mixed 
coupling case — not just the proton-only or neutron-only 
cases — this is indeed the case: the LUX likelihood-based 
limit at 90% CL excludes the entire SD parameter space 
consistent with the DAMA result within the 2a CL, at 
least for the assumed halo model. As always, caveats 
apply. Various assumptions about detector behavior are 
made that, if incorrect, will affect the interpretation of 
the experimental results and alter the WIMP parame¬ 
ter space consistent with those results. See for exam¬ 
ple Ref. [75]. In addition, if the NEST model for low- 
energy events is inaccurate, the low-mass LUX limits may 
weaken. For the more conservative maximum gap analy¬ 
sis, a tiny part of the DAMA-compatible parameter space 
escapes the LUX bounds: On/up = —0.16 ± 0.04 with 
« 10 GeV. This remaining space will be excluded by 
the next LUX data release if excess events are not found. 


B. Application to effective theory 


We now apply the analysis of this paper to a DM model 
for which the SD constraints are particularly important. 
Consider the case of Dirac DM annihilating through the 
s-channel exchange of a spin-1 mediator, V/, via an axial- 
vector interaction. Assuming that the mediator also has 
an axial-vector interaction with SM fermions, we obtain 
the Lagrangian: 


^ D [gxaXl^l^X + 9fafx'"x^f] Vf, 
for DM that is a Dirac Fermion, and 


£ D 


1 

2 


dxaXl^l^X + 


Vn 


(24) 


(25) 


for Majorana DM, where g^a and gfa are unknown cou¬ 
plings [za |5D]. Assuming the limit of low-momentum 
exchange in the WIMP-nucleon scattering process, we 
can integrate out the mediator to obtain the following 
SD scattering cross-section for a Dirac WIMP: 


^{a,a 


irmf, 


J{J+1) 


(Sp) , , (^n) 


J 


Op + 


J 


(26) 
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where J is the spin of the nucleus, is the mediator 
mass and p. is the reduced WIMP-nuclear mass. Here, 
dp and d„ are the WIMP-proton and WIMP-neutron 
couplings, related to the couplings of Section |II A2| via 
= 2\/2GFaN = ^ Majorana 

WIMP, Eqn. (261 adopts an additional factor of 1/2. 
Though we focus here on direct searches for such a par¬ 
ticle, colliders can also place constraints; see Ref. [HD] for 
a discussion. 

Expressions for dp and dn can be derived by start¬ 
ing from the WIMP-quark interactions and performing a 
weighted sum over the components of each nucleon. Since 
the heavy quarks and gluons contribute negligibly to the 
spin content of the nucleon, DM scattering off nucleons 
is dominated by the sum over light quarks. This allows 
us to write: 


OAT = X! 3fa (27) 

q—u,d,s 

where is the nuclear spin content of quark q, and we 
have assumed that the coupling between the WIMP and 
each quark is identical and given by the coupling in our 
above Lagrangian. The standard values for the various 
are A^'^ = A^”^ = 0.84, a 1”^ = A^'’^ = —0.43 and 
A^P’ = Ai"^ = -0.09 0. 

Due to the symmetry in the up and down quark con¬ 
tributions, dn = dp, so On = Up. Neglecting the slight 
difference in the proton and neutron masses, a^p « tr^n 
with 




Xn-^xa 


-1 2 






-\ 2 


E 


_ q—u,d^s 


(28) 


the WIMP-neutron cross-secti(^ obtained using J = 
(5'n) = ^ and (5'p) = 0 in Eqn. (26), and a^p the WIMP- 
proton cross-section. 

For a fixed WIMP mass, a limit on places a limit 
on with that limit corresponding to a linear 


relationship between log{g^a9fa) and logniy, as shown 
in Figure (Figure for a Dirac (Majorana) WIMP. 
Here, we show LUX limits for a 25 GeV WIMP, using the 
90% CL upper limit of a^n < 1-26 x 10“^ pb, determined 
by LUXCalc. We stress a subtle point here: the appro¬ 
priate LUX limit on cr^^n is determined using the Op = On 
relation expected for this model, not the Op = 0 assumed 
in the SD neutron-only case shown in Figure though in 
practice the two limits are similar. We also show in Fig¬ 
ures [^&[^the results of a previous analysis m based on 
the XENONIOO result. As expected, the LUX bounds 


Again, a/v is sometimes used in the literature to refer to the 
normalization used here, as is the case with Ref. ESI- 



FIG. 5. Constraints on the parameters of the effective La¬ 
grangian given in Eqn. pij ) for a 25 GeV Dirac WIMP. The 
LUX SDn limit obtained by LUXCalc is shown in black, with 
coupling values above the black line excluded for a given me¬ 
diator mass. Also shown are the results of a previous analysis 
utilising the XENONIOO SDn limit (green line) |79j The red 
line shows the parameter values required to obtain the cor¬ 
rect DM relic density (Ddm = 0 . 268/0 0 ( 0 ) as measured by 
WMAP and Planck [791181]. 



FIG. 6 . Similar to Figure [^ but for a Majorana WIMP. 


are more stringent, raising the limit on the mass of a 
mediator consistent with relic density observations (red 
curves) from 20 GeV to 30 GeV for a Dirac WIMP, and 
from 10 GeV to 20 GeV for a Majorana WIMP. 


V. CONCLUSIONS 

The particle nature of DM is still unknown. There 
are several theoretical candidates for DM, however the 
WIMP hypothesis remains one of the most popular ex¬ 
planations of the phenomenon. The null results of direct 
detection experiments such as LUX provide limits on the 
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physics of WIMP-nucleus interactions and therefore the 
parameter space of a given WIMP model. 

We have developed and utilized the new tool LUXCalc 
to generate limits on the SI WIMP-nucleon cross-section 
at 90% CL both with Poisson and maximum gap based 
analyses for a WIMP mass in the range [5, 2500] GeV. We 
see that the maximum gap method generally agrees very 
well with the official LUX SI limits, whilst the Poisson 
likelihood-based constraint is weaker by ~ 30% above a 
WIMP mass of about 20 GeV. We then generate the LUX 
SI scattering limits for the low mass region [3,15] GeV 
again using both the maximum gap and Poisson likeli¬ 
hood techniques. For the maximum gap method, we also 
show limits with the progressively more conservative ex¬ 
clusion of all contributions from events with energies be¬ 
low I, 2, and 3 keV. These LUXCalc-generated limits are 
then compared with the official LUX result as well as the 
anomalous signal regions as seen by GoGeNT, CRESST, 
and DAM A and the exclusion curve from SuperCDMS. 
In this mass region, the Poisson-likelihood curve provided 
the strongest constraint. 

We have used LUXCalc to generate for the first time the 
LUX limits on the SD WIMP-nucleon cross-section over 
the full range of WIMP masses. We show constraints 
for both neutron-only (op = 0) and proton-only (on = 0) 
couplings using the maximum gap method detailed in the 
text. We see that for the SD proton-only case, which is 
the SD case where the LUX limits are approximately at 
their weakest, the LUX limit fully excludes the DAM A 
region. In fact, we find that the LUX likelihood-based 
limit at 90% CL excludes the entire SD parameter space 


consistent with the DAMA result within the 2cr CL (for 
the assumed parameters of the SHM). Furthermore, the 
more conservative maximum gap method excludes most 
of the DAMA parameter space, except for the small re¬ 
gion On/op = —0.16 ± 0.04 with « 10 GeV. 

Finally, we have applied the main results of this work 
to an effective theory case where SD constraints on the 
WIMP-nucleon cross-section are particularly important. 
We see that the LUX bounds are more stringent than 
those of a previous study based on the XENONIOO re¬ 
sults m by a factor of ^2 in both the Dirac and Majo- 
rana WIMP cases. 

We have made the LUXCalc tool publicly available for 
future studies involving the LUX results. 
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Appendix A: LUXCalc 


Here we describe the LUXCalc software package, which can be found at Ref. EO] or as ancillary 
files to the arXiv version of this paper. The package provides both a library and a standalone 
program for performing various likelihood and constraint calculations. The software is written in 
Fortran 95, but linking to the library can be easily performed from C-H-. We begin in Section jXT] 
by describing some basic usage of the standalone program, in Section |A 2| we show how to link 
to the library from Fortran, and in Section [A 3| we show how to do the same in C-H-. Finally, in 
Section A 4 we point out a few routines in other software packages (DarkSUSY and micrOMEGAs) 


that may be useful. 

Both the library and program can be compiled by running ‘make all’; however, one of the 
gfortrain or ifort compilers must be installed. 


1. Program 

The LUXCalc program is called in the following form: 

./LUXCalc [mode] [options] [WIMP parameters] 

where [mode] is a flag describing the type of calculation to be performed, [options] are optional 
flags that can be used to set various parameters or control the output, and [WIMP parameters] 
are the WIMP mass and scattering cross-section(s), necessary only in certain modes. Several 
of the program modes are described below and some of the most useful options are described 
in Table jn] A full description of all modes and options can be found by running ‘./LUXCalc 
—help’. 
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—verbosity=l 


—rho=0.4 
—vrot=235.0 


—v0=235.0 
—vesc=550.0 


—Emin=0.0 


Specifies how much detail to provide in the output. Increasing this above the default of 
1 will cause a detailed header to be provided, as well as progressively more detailed data 
output. The additional data provided is mode specific. 

The local density of the DM halo [GeV/cm®]. 

The local rotation speed of the galactic disk [km/s]. The velocity dispersion vo will be 
set to this value as expected for an isothermal spherical halo model (i.e. the Standard 
Halo Model), unless specified via the —vO option. 


The velocity dispersion vg and local escape velocity Uesc (i.e. cutoff speed) [km/s] used 
to define the truncated Maxwellian velocity distribution (Eqn. ( |18| l). Specifically, the 
parameter vg is the most probable speed of the distribution in the absence of any trun¬ 
cation. 

Only consider contributions to the expected signal from events with recoil energies above 
this value [keVj. Detector thresholds are already factored into the efficiencies, so contri¬ 
butions from low-energy events are suppressed regardless of this setting. 


—confidence-level=0.9 


—p-value=0.1 Specifies the confidence level (CL) or the p-value (p = 1 — CL) to use for generating 

constraints. 


—m-tabulation=l.0,1000.0,-20 

For quantities that are tabulated by mass like cross-section limits, this option specifies 
the minimum and maximum masses in the tabulation [GeV], followed by the number of 
tabulation intervals. A negative number for the third value indicates number of intervals 
per decade. 


TABLE 11. Useful options for running the LUXCalc program. Any values shown are the default values. Some options are only 
useful in certain program modes. 


Likelihood. The logarithm of the Poisson-based likelihood, as described in Section Enl is given 
by 

./LUXCalc —log-likelihood [options] [WIMP parameters] 

The WIMP parameters are a list of values: the WIMP mass [GeV] followed by one, two, or four 
cross-sections [pb]. In the first case, the single cross-section is the SI WIMP-nucleon cross-section, 
assumed to be the same for protons and neutrons. In the second case, the two cross-sections 
are the SI and SD WIMP-nucleon cross-sections, again assuming identical couplings for protons 
and neutrons. In the last case, the four cross-sections are, in order, the SI WIMP-proton, SI 
WIMP-neutron, SD WIMP-proton, and SD WIMP-neutron cross-sections. 


Maximum gap p-value. The maximum gap p-value, as described in Section m is given by 


./LUXCalc —log-pvalue [options] [WIMP parameters] 


Specifically, this returns the quantity 1 — Co, where Co is defined in Eqn. (23) and is technically 
only an upper limit on the p-value, not the p-value itself. As opposed to the Poisson-based 
likelihood above, this analysis involves no background subtraction and is hence conservative. 
The WIMP parameters are as described for the likelihood mode above. 


Likelihood constraints. Tabulated (by mass) upper and lower cross-section constraints (i.e. 
confidence intervals) as determined via the likelihood are generated by 

./LUXCalc —constraints-Sl [—theta-Sl-pi=0.25] [options] 

./LUXCalc —constraints-SD [—theta-SD-pi=0.25] [options] 

where the two lines correspond to spin-independent (SI) and spin-dependent (SD) interactions. 
When determining constraints, the ratio between WIMP-neutron and WIMP-proton couplings 
will be kept fixed, with the ratio defined in terms of the polar angle 0 in the Gn-Gp plane, 
i.e. tan0 = GnjGp = /n//p = anjap. The default behavior is to take the two WIMP-nucleon 
couplings to be identical; otherwise, 0 can be specified (in units of tt) via the options as shown. In¬ 
crease the verbosity (e.g. — verbosity=3) to show the corresponding WIMP-neutron constraints 
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in addition to the WIMP-proton constraints. The confidence level (CL) of the confidence inter¬ 
vals is specified via the —confidence-level option; the default is 90% CL. 

Maximum gap limits. Tabulated (by mass) upper limits on the cross-section(s) as determined 
by the maximum gap method are generated by 

./LUXCalc —limits-Sl [—theta-Sl-pi=0.25] [options] 

./LUXCalc —limits-SD [—theta-SD-pi=0.25] [options] 

where the options are as described for the likelihood constraints above. 


2. Library: Fortran usage 

LUXCalc is written as a single, self-contained Fortran 95 module. All floating point values are 
in the REALMS format, while integers are of type INTEGER. The module must be loaded in any 
user routine that calls LUXCalc routines: 

USE LUXCalc 


Initialization. Before calling any routines, the module must first be initialized with 
CALL LUXCalc_lnit(intervals=.TRUE.) 

The single argument here specihes if calculations should be performed for the intervals (gaps) 
between events. This is necessary for generating maximum gap limits, but is not required for any 
likelihood calculations. This initialization need only be performed once. To force only recoils of 
energy greater than ifmin to be considered in calculating rates, use 

CALL LUXCalc_SetEmin(Emin=OdO) 

where the argument is in keV. As noted elsewhere, detector thresholds are already factored into 
the efficiencies, so contributions from low-energy events are suppressed regardless of this setting. 

Halo model. The parameters of the Standard Halo Model can be specified via 

CALL LUXCalc_SetSHM(rho=0.4d0,vrot=235d0,v0=235d0, & 
vesc=550d0) 


where the arguments are the local DM density po [GeV/cm^], the disk rotation speed Urot [km/s], 
the velocity dispersion vq [km/s], and the galactic escape speed Vesc [km/sj. The values shown 
here are the defaults, which are already set when LUXCalc_Init () is called, so the above function 
call is spurious. The halo parameters can be modified at any time. 


WIMP mass and couplings. The WIMP parameters are specified by one of three routines: 

CALL LUXCalc_SetWIMP_mfa(m,fp,fn,ap, an) 

CALL LUXCalc_SetWIMP_mG(m,GpSI,GnSI,GpSD,GnSD) 

CALL LUXCalc_SetWIMP_msigma(m,sigmapSI,sigmanSI, & 

sigmapSD,sigmanSD) 


Here, m is the WIMP mass [GeVj; fp and fn are SI WIMP-proton and WIMP-neutron 
couplings [GeV“^], respectively; ap and an are SD WIMP-proton and WIMP-neutron couplings; 
GpSI, GnSI, GpSD, and GnSD are WIMP-nucleon couplings [GeV“^], differing from / and a only in 
normalization as discussed in Section HA and the sigma arguments are WIMP-nucleon cross- 
sections [pbj. A negative value for a cross-section can be used to indicate the corresponding 
coupling should be taken to be negative. The current mass, couplings, and cross-sections can be 
retrieved with the corresponding routines 


CALL LUXCaIc_GetWIMP_mfa(m,fp,fn,ap,an) 

CALL LUXCaIc_GetWIMP_mG(m,GpSI,GnSI,GpSD,GnSD) 

CALL LUXCaIc_GetWIMP_msigma(m,sigmapSI,sigmanSI, & 

sigmapSD,sigmanSD) 
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The returned cross-section values will not be set negative for negative couplings. 

Calculations. After any changes to the WIMP parameters and/or the halo distribution, the 
LUX rate calculations must be performed using 

CALL LUXCalc_CalcRates0 


This routine performs the various LUX rate calculations that are used for determining likelihoods 
and the various constraints. Thus, this routine must be called before obtaining expected events, 
likelihoods, p-values, etc. The relevant quantities are stored internally. 

Events. The number of observed and expected events for the current WIMP are provided by 
the functions: 

INTEGER FUNCTION LUXCaIc_Events() 

REAL*8 FUNCTION LUXCaIc_Background() 

REAL*8 FUNCTION LUXCaIc_SignaI() 

REAL*8 FUNCTION LUXCaIc_SignaISI0 
REAL*8 FUNCTION LUXCaIc_SignaISD() 

In order, these return the observed number of events in LUX, the average expected background 
events, the average expected signal events, and the separate SI and SD contributions to the 
expected signal. 


Likelihoods and p-values. The statistical functions for evaluating LUX results in the context 
of the current WIMP are: 

REAL*8 FUNCTION LUXCaIc_LogLikeIihood() 

REAL*8 FUNCTION LUXCaIc_LogPVaIue() 

REAL*8 FUNCTION LUXCaIc_ScaIeToPVaIue(logp) 


The first function returns the log of the likelihood using a Poisson distribution in the number 
of observed events given the expected background and signal. The second function returns 
the logarithm of the p-value in a more conservative no-background-subtraction analysis. If 
LUXCaIc_Init0 was called with a .TRUE, argument, then the maximum gap method is used, 
with p = 1 — Co, where Co is given by Eqn. (23). Otherwise, A Poisson distribution with zero 
background contribution is assumed. This function is only intended for determining conservative 
one-sided limits as the value returned is technically an upper limit on the p-value and not the 
p-value itself. The last function determines such a limit by identifying the factor x such that 
(T = xa-Q gives the desired p-value (specified as log(p)), with erp the currently specified WIMP- 
nucleon cross-sections. The quantity xao is then the limit on the cross-sections at the given 
confidence level (1 — p), assuming a fixed ratio of WIMP-nucleon couplings (e.g. /p = /„). 


3. Library: C-|—h usage 

To make usage of LUXCalc with C-h- code easier, a C-h- interface file LUXCalc .hpp is provided. 
The routines and functions are the same as those described for Fortran in the previous section, 
with identical names and signatures, though Fortran subroutines become C-h- void functions. 
All arguments and return values are bool, int or double. 


4. Useful software 

Here we show how to extract the relevant DM parameters from two of the most popular 
software packages for examining DM in the context of SUSY: DarkSUSY [T71[5B] and micrOMEGAs 
[T8l [57U5^ . DarkSUSY is written in Fortran 77, with the various G WIMP-nucleon couplings for 
a given SUSY model provided by the dsddgpgn routine. The WIMP (neutralino) mass must 
be retrieved from various common blocks. The necessary parameters can be retrieved from 
DarkSUSY and set in LUXCalc as shown in this minimal Fortran 95 example: 
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/ Load LUXCalc module 
USE LUXCALC 

! Variables to store WIMP mass and couplings 
REALMS :: M,GpSI,GnSI,GpSD,GnSD 
/ Calculated quantities (examples) 

REAL*8 :: Inlike,signal 

/ DarkSUSY common blocks defined in ’dsmssm.h’ 

CHARACTER*8 :: pacodes_ctmp(0:50) 

INTEGER :: pacodes_itmp(60),kn(4),lsp,kln 
REAL*8 :: mass(0;50),mspctm_rtmp(6) 

COMMON /PACODES/ pacodes_itmp(l:18),kn,pacodes_itmp(23:60), & 
pacodes_ctmp 

COMMON /MSPCTM/ mass,mspctm_rtmp 
COMMON /MSSMIUSEFUL/ Isp.kln 

/ Initialize LUXCalc 
CALL LUXCalc_lnit(.FALSE.) 

! For each SUSY model, do following »»»» 

! Set WIMP parameters 
M = mass(kn(kin)) 

CALL dsddgpgn(GpS1,GnS1,GpSD,GnSD) 

CALL LUXCalc_SetWlMP_mG(M,GpSl,GnSl,GpSD,GnSD) 

! Calculate rates for current WIMP 
CALL LUXCalc_CalcRates() 

! Get likelihood, expected signal events, etc. 

Inlike = LUXCalc_LogLikelihood() 
signal = LUXCalc_Signal() 


The micrOMEGAs package provides the WIMP mass and couplings in the Mcdm global variable 
and nucleonAmplitudes routine, respectively; see Ref. |18j for a description of the relevant 
micrOMEGAs coupling routine and its arguments. A minimal C++ example using micrOMEGAs is: 

// Initialize LUXCalc 
LUXCalc_lnit(false); 

// For each SUSY model, do following »»»» 

// Set WIMP parameters 

// separate particle/anti-particle couplings 
double lambdap[2],lambdan[2],xip[2],xin[2]; 

// FeScLoop is micrOMEGAs-provided function 

nucleonAmplitudes(FeScLoop,lambdap,xip,lambdan,xin); 

double M = Mcdm; // Mcdm is global variable 

double GpSl = 2*lambdap[0]; 

double GnSl = 2*lambdan[0]; 

double GpSD = 2*xip [0]; 

double GnSD = 2*xin[0]; 

LUXCalc_SetWlMP_mG(M,GpSl,GnSl,GpSD,GnSD); 

// Calculate rates for current WIMP 
LUXCalc_CalcRates(); 

// Get likelihood, expected signal events, etc. 
double Inlike = LUXCalc_LogLikelihood(); 
double signal = LUXCalc_Signal(); 
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We finally point out one subtlety: DarkSUSY and micrOMEGAs will yield somewhat different 
WIMP-nucleon scattering cross-sections for a given SUSY model due to different values chosen 
for the hadronic matrix elements that go into calculating WIMP-nucleon couplings from the 
WIMP-quark couplings [HI]- This is more of an issue for the SI cross-sections, which typically 
differ by a factor of 0{2). Both packages allow the matrix elements to be modified; see their 
respective manuals for further details. 
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